Probing coherent charmonium photoproduction off light nuclei at medium energies 
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We demonstrate how the elementary amplitudes jN — > , the amplitude of the nondiagonal 
J/ipN ip'N transition, and the total J/ipN and ip'N cross sections can be determined from 
measurements of the coherent J/ip and ip' photoproduction off light nuclei at moderate energies. For 
C**) ■ this purpose we provide a detailed numerical analysis of the coherent charmonium photoproduction 

off silicon within the generalized vector dominance model (GVDM) adjusted to account for the 
physics of charmonium models and color transparency phenomenon. 

I. INTRODUCTION 

First measurements of quasielastic photoproduction of J/ip and ip' mesons off nucleon and nucleus targets at 
moderate energies were performed in the seventies at SLAC and Cornell [1-3] . The measurement of the A-dependence 
of J/ip photoproduction [2] was used by the authors to extract the J/ip - nucleon absorption cross section a(J/ipN) — 
3.5 ±0.8 mb, using a semi-classical model for J/ip - nucleus interactions. At the same time the estimates of a(J/ipN) 
C^ - • based on the Vector Dominance Model and the optical theorem yield a(J/ipN) = 1.3 ± 0.3 mb and a(ip'N) w 
0.7a(J/iPN) [1]. This small value of the J/ipN cross section comparing to that for irN interaction is in variance with 
the preQCD ideas on universal hadron-hadron interaction. On the contrary, these observations can be considered [4] 
as evidence for the previously predicted [5,6] genuine QCD effect of the dependence of the cross section on the size 
of the region occupied by color within a hadron. Besides, QCD predicts that a^'N should be significantly larger 
■ then o~j/^,n because the radius of the ip' is two times larger than the radius of the J/ip. Recently first data on the 
hadroproduction of J /ip and ip' in pA collisions confirmed these QCD predictions but large experimental errors and 
the complicated theoretical analysis preclude an unambiguous extraction of the charmonium-nucleon cross sections 
from the data. As a result the general situation with the extraction of these cross sections is still unsatisfactory - the 
Oh, determined values of o~j/^ N range from 2 mb to 6 mb and, correspondingly, the range of a^N is from 8 mb to almost 
20 mb [7]. On the other hand the precise values of these cross sections are urgently needed for the unambiguous 
interpretation of the yield of particles with hidden charm produced in the ultrarelativistic heavy ion collisions (for 
the review and references see e.g. [8]). 

Since direct measurements of the charmonium-nucleon cross sections are impossible, the photoproduction of the 
J/ip and the ip' off nuclei at coherent moderate energies appears to be one of the most promising tool to measure these 
quantities. Really, in the moderate energy domain both the coherence length and the formation length are smaller or 
comparable to the radii of nuclei. So, a photon produces a cc pair which transforms to a charmonium state inside the 
nucleus with a noticeable probability and this meson can interact with the surrounding nucleons. However, it appears 
that the standard analysis of the data within the Glauber based Vector Dominance Model (VDM) which provides a 
more or less reasonable treatment of the light vector meson photoproduction at low and moderate photon energies 
fails in the case of the charmonium photoproduction. It was understood that according to the charmonium models 
the very small size of the cc pair which is controlled in the production point by the mass of the c-quark leads to the 
radical change of the predictions based on the VDM [9-11]. 

The aim of this paper is to suggest a new option for the extraction of charmonium-nucleon amplitudes from 
measurements of the coherent charmonium photoproduction off light nuclei at moderate photon energies. In a recent 
paper [12] we demonstrated that large masses of the charmonium states (my) and the fluctuations of the strength 
of the interaction of different charmonium states with nucleons inside the nuclear medium expected in QCD lead to 
significant cross section oscillations in this process. We are interested in the photon energy range 15 GeV < E 1 < 

2 

40 GeV relevant for the planned SLAC experiment [13]. In this case the minimal momenta w transferred 
to nucleus are 0.35 GeV > q^ > 0.15 GeV and the significant oscillations due to the oscillating nuclear form factor 
are spectacularly revealed [12] in the energy dependence of the coherent charmonium photoproduction cross sections. 
Another source of oscillations predicted by the QCD color screening phenomenon is the significant difference of the 
amplitudes of J/ipN — > J/ipN, ip'N — ► ip'N, and J/ipN <-> ip'N which leads to oscillations of the rescattering strengths 
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and noticeably modifies the energy dependence of the cross sections. We suggest to explore these oscillations as a new 
method of determining the elementary photoproduction amplitudes and the amplitudes of J/ip and ip' interactions 
with nucleons. 

Experimentally the suggested approach can be realized at SLAC where measuring of the charmonium photopro- 
duction at low energies (15 - 35 GeV) is planned [13]. 



II. DESCRIPTION OF THE MODEL 



The key feature of the charmonium photoproduction is that the cc configuration of the photon wave function is 
spatially small in the production point because its transverse size is controlled by the large mass of the c-quark. 
However, at the photon energies of 15 -j- 40 GeV wc arc interested in, the charmonium formation length, If « 
uj[2m c (m^,i — mj^)] -1 , is still comparable to the internucleon distance and the nucleus radius. Hence, the noticeable 
probability exists that at these moderate energies a spatially small configuration would transform to hadron states 
before a collision with a second nucleon. In this case using a hadronic basis which properly takes into account 
nonperturbative QCD effects seems to be relevant for the calculation of the photoproduction cross section 
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Here \I r = J/ip rcsp. ip , —t m i n = is the longitudinal momentum transfer in the 7 — > ^ transition. The 
photoproduction amplitude F~,A-*VA(t) is normalized so that ct 7 a->*a(w) is in lib. To calculate this amplitude the 
multistcp production Glauber approach formulae [14] were combined with the GVDM [15] which we adjusted for the 
charmonium case in order to for account the color screening effects. A more detailed description of this Generalized 
Glauber Model (GGM) (which we successfully tested [16] in the case of the p-meson coherent photoproduction) is given 
in ref. [12]. Here we briefly summarize several essential points. Since the coherent production of heavy charmonium 
states off nuclei at low and moderate energies is suppressed by the target form factor it is legitimate to build the GVD 
model restricted to two states - J/ip and ip': 

f-fN^J/i/jN = 7 fj/ipN^J/ipN + -7—ftj)'N^J/ij)N j (2) 

Jj/tp H' 

6 6 

fjN^ip'N — -J—fip'N^ip'N + 7 fj/iPN^ip'N ■ (3) 

Jif>' J J/ip 

The charmonium-nucleon coupling constants — 10.5 ± 0.7, and — 30.9 ± 2.6 are determined from the widths 
of the vector meson decays * — ► ee. Since in the photoproduction processes cc pair is produced within spatially 
small configuration one can neglect for a moment the direct photoproduction amplitude and obtain from eqs. (2,3) 
the approximative relations between amplitudes 

fiP'N^tp'N ~ — fiP'N^J/ipN ~ f t fj/rpN^J/rpN ■ (4) 

JJ/ip J J/ip 

Note that a comparatively large absolute value of the nondiagonal amplitude and the negative value of the ratio of 
nondiagonal and diagonal amplitudes follow from the presence of the color screening in the charmonium photopro- 
duction as dictated by QCD. For a more accurate determination of the elementary ip'N —> ip'N and J/ipN <-> ip'N 
amplitudes within GVDM (for the details see [12]) we use as an input in eqs. (2,3) the absolute values of the forward 
7A — > J/xpN cross section, the relation da 7 N^ip' N{t m in) /dt = 0.15-dcr 7 Ar^j/^,7v(tmin)/di and the J/ipN cross section 
cj/ipN ~ (3-5 ± 0.8) mb found at moderate energies at SLAC [1,2]. To show explicitly how uncertainty in the value 
°f vj/ipN influences on the photoproduction cross section we used estimates of amplitudes determined within the 
corridor of the experimental errors. This is shown by the filled area in all results of calculations. 

We performed calculations within the Generalized Glauber Model (GGM). Note however, that in the case of coherent 
photoproduction of charmonium off light nuclei at moderate energies, the parameter pa^^RA is small and one can 
simplify calculations by taking into account only one rescattering of the produced heavy meson. Then only two 
forward amplitudes provide the dominating contribution and the coherent charmonium photoproduction is described 
by the sum of the diagrams in fig. 1. The first graph is the direct photoproduction amplitude while the second one 
accounts for one rescattering of the produced charmonium states within the nuclear medium. 
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FIG. 1. Diagrams describing the coherent charmonium photoproduction: left - direct photo-production amplitude, right - 
photoproduction with one rescattering . 

Since V,V = J/ip,ip' the second diagram includes the diagonal as well as nondiagonal rescatterings. We want to 
emphasize that in the case of nondiagonal transition the effect of charmonium absorption enters in the next to leading 
order. 



Within the Generalized Glauber Model we performed calculations for a wide range of nuclei trying to choose the 
optimal nuclei for determining the elementary amplitudes from measurements of the coherent photoproduction of the 
J/ip and the tp'. The rates are small for the lightest nuclei while for sufficiently heavy nuclei multiple rescattering 
effects become noticeable. Hence we found that nuclei with A = 20 -f- 40 are optimal for our purposes. Thus in the 
following we will present the results for silicon nuclei. The silicon target is also of interest since one could envision an 
active silicon target which would allow to select the coherent events in a cleaner way. 

The energy dependence of the J/ip cross section is compared (fig. 2) to that obtained within the Impulse Approxi- 
mation where 



Skyrme (HFS) model which provided a very good (with an accuracy w 2%) description of the global nuclear properties 
of spherical nuclei along the periodical table from carbon to uranium [17] and the shell momentum distributions in 
the high energy (p,2p) [18] and (e,e'p) [19] reactions. 



III. RESULTS AND DISCUSSION 
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FIG. 2. Energy dependence of the forward coherent 7 + Si — > J/i/> + S'i photoproduction cross section calculated in the 
Generalized Glauber Model compared to the cross sections in the Impulse Approximation (dotted line). The forward elementary 
7./V — > J/ipN cross section is shown by the dashed line. 
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FIG. 3. Energy dependence of the ratio of the forward coherent 7 + Si — > J/tp + Si photoproduction cross section calculated 
in the Generalized Glauber Model to that in the Impulse Approximation. 
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The distinctive feature of the coherent charmonium photoproduction cross sections in the kinematics we are inter- 
ested in is that the dependence of the cross section on the photon energy is governed by the longitudinal nuclear form 
factor at a relatively large value of t min which becomes more important with a decrease of the photon energy. We 
want to emphasize that the use of a realistic nuclear density in our calculations ensures a reasonable description of 
the nuclear form factor in the relevant momentum transfer range. One can infer from fig. 3 that the cross section of 
the J/tp production off the silicon calculated in GGM practically coincides with the Impulse Approximation result. 
The deviation does not exceed 4 5% and it varies weakly when the rescattering amplitudes are changed within 
region allowed by the uncertainties of the input values of the J / ipN cross section as determined at SLAC [2] . This 
implies that one should have an unprecedented accuracy both in the measurements and in the theoretical analysis 
in order to extract the J/ipN elementary cross section from such measurements. However these measurements can 
be used to determine precisely the elementary forward photoproduction cross section. Really, the contribution of the 
diagonal J/i/jN rescattering can be neglected. Then the comparison of the data with the cross section calculated 
in the Impulse Approximation (eq. 5) with the nuclear form factor which is well determined from the high energy 
electron- nucleus elastic scattering immediately provides us with the forward elementary — > J/ipN cross section. 
As it is seen from fig. 2, the low t cross section of J/tp production off the silicon in spite of the strong suppression by 
the nuclear form factor is considerably larger than the der^jy— >j7i/;iv {t m in)/dt (dashed line in fig. 2). Hence, one may 
significantly improve the accuracy of the determination of d<7 7 jv^j/^./v(imm)/di. 
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FIG. 4. The energy dependence of the forward coherent 7 + Si — » tp' + Si photoproduction cross section calculated in the 
Generalized Glauber Model (solid lines, filled blue) compared to the cross section in the Impulse Approximation (dotted line) 
and to the cross section calculated in GGM without accounting for the diagonal ipN — > ipN rescattering (dashed lines, yellow). 
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FIG. 5. The energy dependence of the ratio of forward coherent ip'-to-J /ip photoproduction cross sections calculated in the 
Generalized Glauber Model compared to the ratio in the Impulse Approximation (dotted line), to the ratio of cross sections 
calculated in GGM without accounting for the diagonal rescatterings (dashed line, yellow filled) and to the ratio without 
accounting for the direct ip' photoproduction (dashed lines, violet filled) . 

The situation is qualitatively different for the coherent ip' photoproduction off silicon (fig. 4). The cross section 
calculated within GGM significantly exceeds the Impulse Approximation cross section. This is because the suppression 
due to the diagonal rescatterings is negligible as compared to the enhancement due to the contribution of the two-step 
nondiagonal — -> J/tpA — > ip' A process. The range of energy in the vicinity of the form factor minimum for the 
forward ip' production is especially interesting. In this case the main contribution to the cross section (fig. 5) originates 
from the nondiagonal rescattering. In particular, at u> w one can extract from the data the nondiagonal 

elementary J/ipN — > ip'N amplitude by measuring the forward relative ip' -to- J/ip yield. In such a ratio all other 
inputs, namely, the nuclear density and the elementary 7 AT — > J/ipN amplitude are already fixed. Moreover, since 
they enter both in the numerator and denominator, the major uncertainties are canceled out. 

In the discussed energy range the soft rescattering processes within the nuclear medium are characterized 
by a rather weak energy dependence, like s 08 , which can be neglected to a very good accuracy in the ratio of the 
forward ip' and J/ip photoproduction cross sections. Hence the energy dependence of this ratio (fig. 5) should originate 
primarily due to the contribution of the direct ip' production. Thus one would be able to determine the 7 A — > tp'N 
amplitude by measuring the relative ip'-to-J/ip yield. 

We want to emphasize that the suggested procedure for extracting the nondiagonal amplitude and amplitudes of 
direct J/ip and ip' photoproduction from the nuclear measurements is practically model independent. 

The task of determining the diagonal J/ipN —* J/ipN and tp'N — > ip'N amplitudes which are relevant for the 
suppression of the charmonium yield in heavy ion collisions is much more complicated. 

If one would naively use the VDM, neglecting nondiagonal transitions and using values of a to t{ J/ipN) based on the 
SLAC data [2] one would predict a rather significant suppression of the J/ip yield: » 10% -j- 15% for light nuclei, and 

30% -^40% for heavy nuclei. However, we find [12] a strong compensation of the suppression due to the contribution 
of the nondiagonal transitions. As a result we find that overall the suppression does not exceed 5% -^10% for all nuclei 
along the periodical table. Hence, an extraction of the diagonal amplitudes from the measured cross sections would 
require a comparison of high precision data with very accurate theoretical calculations including the nondiagonal 
transitions. Thus one may try another strategy. If the elementary -fN — > J/ipN, and -fN — > ip'N photoproduction 
amplitudes, as well as the nondiagonal amplitude J/ipN <-> ip'N would be reliably determined from moderate energy 
data it would be possible to determine the imaginary parts of the forward diagonal amplitudes from the GVDM 
equations (2,3) reversing the procedure which we used in section 2, namely 
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(6) 



Qfy'N^Hp'N — 1 




^fj/ipN->ip>N 



(7) 



The main limitation of the suggested procedure is the restriction of the hadronic basis to the two lowest IS, 2S 
charmonium states with the photon quantum numbers: J/ip , tp'. This question was discussed in [12]. The disregard 
by the closest higher charmonium state ip(3770) can change the estimate of the J/tp N amplitude by « 10%. Influence 
of the higher mass resonances is expected to be even weaker - the constants 1/fv relevant for the transition of a 
photon to a charmonium state V rapidly decrease with the resonance mass. This is because the radius of a bound 
state, ry, is increasing with the mass of the resonance and therefore the probability of the small size configuration 
being oc l/r v is decreasing with an increase of mass (for fixed S,L). Besides, the asymptotical freedom in QCD 
dictates decreasing of the coupling constant relevant for the behavior of the charmonium wave function at small 
relative distances. Experimentally one finds from the data on the leptonic decay widths that 1/fv drops by a large 
factor with increasing mass. An additional suppression arises due to the weakening of the soft exclusive nondiagonal 
VN <-> V'N amplitudes between states with the different number of nodes. Hence, determining the imaginary parts 
of diagonal rescattering amplitudes from the GVDM equations seems to be possible. Since in the medium energy 
domain the energy dependence of soft rescattering amplitudes is well reproduced by a factor s 08 the real parts can 
be found using the well known Gribov-Migdal relation 



Using the suggested procedure one will be able to determine all elementary amplitudes with a reasonable precision 
by measuring the photoproduction of J/tp and tp' off the light nucleus at the medium photon energies. The cross 
check of this approach would be a comparison of the cross sections calculated within GGM with parameters fixed in 
the analysis of the light nuclei with the experimental cross sections measured in photoproduction of charmonium off 
heavy nuclei as well as in the quasielastic processes. 



We used the Generalized Glauber Model which combines the multistep Glauber approach and the Generalized 
Vector Dominance Model to calculate the coherent charmonia photoproduction off light nuclei at medium energies 
corresponding to the kinematics of the planned E160 SLAC experiment. We demonstrate that the nondiagonal 
amplitudes which follow from the significant QCD color fluctuations within hadrons give a significant contribution to 
the coherent cross section of the tp' photoproduction off light nuclei in this energy domain via the two step production 
mechanism, while the production of the J/tp is well described by the Impulse Approximation. We show how one 
can determine the elementary amplitudes of the charmonia photoproduction as well as the genuine J/ipN, tp'N cross 
sections and the nondiagonal amplitude from measurements in this domain. 
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